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Abstract 

Rotating accretion flow may undergo centrifugal pressure mediated shock transi- 
tion even in presence of various dissipative processes, sucli as viscosity and cooling 
mechanism. The extra thermal gradient force along the vertical direction in the 
post shock flow drives a part of the accreting matter as bipolar outflows which are 
believed to be the precursor of relativistic jets. We compute mass loss rates from a 
viscous accretion disc in presence of synchrotron cooling in terms of the inflow pa- 
rameters. We show cooling significantly affects the mass outflow rate, to the extent 
that, jets may be generated from flows with higher viscosity. We discuss that our 
formalism may be employed to explain observed jet power for a couple of black hole 
candidates. We also indicate that using our formalism, it is possible to connect the 
spectral properties of the disc with the rate of mass loss. 
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1 Introduction 



In recent years, it has been established that AGNs and Microquasars suffer 
mass loss in the form of jets and outflows (IFerraril . Il998l : iMirabel fc Rodriguezl . 
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19991 ). Generation of jets or outflows around gravitating centres with hard 
boundaries {e.g., neutron stars, YSOs etc.) are quite natural, however, it is 
altogether a different proposition to consider the same around a black hole. 
As black holes do not have either hard boundaries or intrinsic atmospheres, 
jets/outflows have to originate from the accreting matter onto black holes, 
though there is no consensus about the exact mechanism of jet formation. One 
of the motivation of studying black hole accretion is therefore to understand 
the primary mechanism in the accretion process which may be responsible 
for the generation of jets. In addition, recent observations have established 
that, whatever be the exact mechanism behind the formation of jets/outfiows 
around black holes, the formation of jets is intrinsically linke d with spec 



tral st ates of the associated black hole candidates. In particular, iGallo et al. 



(I2OO3I ) showed that quasi steady jets are generally ejected in the hard state, 
which suggests that the generation or quenching of jets do depend on various 
states of the accretion disc. Several theoretical attempts were ma de to explain 
the p ossible mechanisms of jet generation from accretion disc. IXu fc Chen 
(119971) reported the formation of outflows by cons idering self-similar solu- 



tions. 



Chakrabartil (119991 ): iDas fc Chakrabartil (119991 ) estimated mass outflow 
rates in terms of inflow parameters from an inviscid advective disc. In par- 
ticular, these authors showed that the centrifugal barrier may produce shock, 
and the post-shock disc can generate bipolar outflows. They also showed mass 
outflow rates depe nd on the s t rength of the centrifugal barrier, as well as, 
its thermal driving. iDas et al. 1 (l2001bl ) extended this work to show that such 
outflows generated by accretion shock is compatible with the spectral state of 
the accretion disc. The shock induced relativistic outflows could be obtained if 
variou s acceleration mechanism, namely, first ord er Fermi acceleration a t the 
shock (ILe fc Beckerl . |2005| ). or radiation pressure (IChattopadhyayl . l2005l ). are 
considered. 



Recently, IChattopadhyay &: Dad (120071 ) computed mass outflow rates from a 
viscous advective disc and showed that the mass outflow rate decreases with 
the increase of viscosity parameter. In realistic accretion disc, a variety of 
dissipative processes are expe cted to be present, and viscosity is just one of 



TO 

them. In absence of mass loss, iGu fc Lul (120041 ) conjectured that cooling pro- 
cesses will not affect the nature of advective accretion solutions. However, 



Dasi (I2OO7I ) explicitly showed that cooling processes play a crucial role in de- 
termining the flow variables as well as the shock properties. Therefore, it will 
be worthwhile to investigate, how cooling would affect the mass outflow rate 
from a viscous accretion disc. In presence of viscosity, as matter flows inward 
angular momentum decreases while specific energy increases. A cooling pro- 
cess unlike viscosity, only reduces the energy of the flow and leaves the angular 
momentum distribution un-affected. Thus the increase of flow energy due to 
viscous heating may be abated by incorporating cooling mechanism. As cool- 
ing is more efficient at the hotter and denser post-shock region (abbreviated as 
CENBOL = CENtrifugal pressure supported BOundary Layer), the decrease 
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of CENBOL energy will be more pronounced compared to the pre-shock en- 
ergy. In reality, more energetic flows at the outer edge, which do not satisfy 
shock conditions in absence of cooling, may undergo shock transition in its 
presence. Consequently, more energetic CENBOL may be produced for flows 
with higher cooling efficiency, and hence there is a possibility of enhanced jet 
driving. In this paper, we would like to address these issues in detail. 

In the next section, we present the model assumptions and the governing equa- 
tions. In Section 3, we discuss the methodology of computing self-consistent 
inflow-outflow solutions and present the solutions. In Section 4, we apply our 
formalism on two black hole candidates to compute the mass outflow rate, and 
compare it with the observed jet power. In the last section we draw concluding 
remarks. 



2 Model Assumptions and Equations of motion 



In a disc-jet system, there are two separate flow geometries, namely, one for 
accretion flows and the other for outflows. Axis-symmetry and steady state 
conditions are assumed for the disc-jet system. In the present paper, we con- 
sider thin, viscous accretion flow in presence of synchrotron cooling. Jets are 
assumed to be tenuous. Since jets are in general collimated, they should have 
less angular momentum and therefore less differential rotation compared to 
the accretion disc. Thus, we ignore the effect of viscosity in jets. As jets are 
believed to originate from the inner part of the disc, which in our model is 
the CENBOL, the jet base must be described by identical local accretion flow 
variables (see section 3), i.e., the specific energy, the angular momentum etc 
of the CENBOL. Consequently, we neglect the torque between the disc and 
the jet at the jet base. It is to be remembered that, to keep the jets collimated, 
angular momentum will be reduced either by magnetic field (stochastic fields, 
considered in th e pape r, are not effective in doing so), or by radiation [see. 



(IChattopadhyayl . l2005l )]. however these processes have not been considered 
here. In reality, back reactions on the disc in the form of extra torque at the 
jet base and/or feedback effect from failed jets are not altogether ruled out. To 
study these effects, one requires to undertake numerical simulation, which is 
beyond the scope of the present frame work. Moreover, jets are supposed to be 
colder than the accretion discs. Therefore, we assume jets to be adiabatic, at 
lea st up to its critical point. We use pseudo- Newtonian potential introduced 



by iPaczyhski fc Wiital (Il980l ) to approximate the space time geometry around 



a non-rotating black hole. 

A schematic structure of shocked advective accretion disc and the associated 
jet are presented in Fig. 1. Here, Xco and the outer and the inner critical 

points of the disc, respectively. The centrifugal pressure acts as a 'barrier' to 
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Fig. 1. A schematic diagram of disc-jet system. The outer and inner critical points 
Xco and Xci are marked in the figure. The shock is located at Xg- The jet geometry 
is bounded by FW and CB. MM' = xpw and MM" = xcb (described in the text). 

the supersonic matter at Xd < x < Xco and a shock at Xg is formed. The post- 
shock disc is indicated in the figure as CENBOL. At the shock, matter mo- 
mentarily slows down and ultimately dives into the black hole supersonically 
through Xci- Excess thermal driving in CENBOL drives a fraction of accret- 
ing matter as bipolar jet which flows within two geomet ric surfaces called the 



Funne l Wall (FW) and the Centrifugal Barrier (CB) (IMolteni et al. 
1996al ). 



1994 



The system of units used in this paper is 2G = Mbh = c = 1, where G, Mbh 
and c are the universal gravitational constant, the mass of the black hole and 
the speed of light, respectively. Since we use the geometrical system of units, 
our formalism is applicable for both the galactic and the extra galactic black 
hole candidates. Two separate sets of hydrodynamic equations for accretion 
and jet, are presented bellow. 

The dimensionl ess hydrodynamic equations th at govern the motion of accret- 
ing matter are (IChakrabartil . Il996l : iDad . 120071 1 , 



the radial momentum equation : 



dx p dx x'^ 



2(x 



(la) 



where, u, p, P, and A(x) are the radial flow velocity, the local density, the 
isotropic pressure and the local specific angular momentum, respectively. Here 
X is the cylindrical radial coordinate. 
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The baryon number conservation equation 



(16) 



where, M and S are the mass accretion rate and the vertically integrated 
density, respectively. In our model, the accretion rates in the pre shock and 
post shock regions are different as some fraction of the accreting matter is 
ejected as outflow. Actually, the post-shock matter is flown into two channels 
— one is the accreting p art (falling onto black holes through x^) and the other 



is th e outflowing part (IMolteni et al. I . Il994j . Il996al : IChattopadhyay &: Das 



20071 ). More specifically, the combination of accretion and outflow rate in the 
post shock region remain conserved with the pre-shock accretion rate (see Eq. 
3). 



The angular momentum conservation equation : 



u- 



dx 



1 d f 



TjX dx ^ 



0, 



(Ic) 



where, Wx(i,{= —all) denotes the viscous stress, a is the viscosity parameter 
and n is the vertically integrated total {i.e., thermal+ram ) pressure. The vis- 
cosity prescription employed in this paper was developed by lChakrabarti &: Molteni 
( I1995I ) a nd has been e r nploy e d to study advective accre t ion disc by a group of 



worke rs dChakrabartil. Il996l: IChakrabarti fc Da.4 l2004l : IGu fc Lil booi [Pas 



2OO7I : IChattopadhvav fc Dad . l2007h . This viscosity prescription is more suit- 
able for flows with significant radial velocity as it maintains angular momen- 
tum distribution continuous across the shock unlike Sakura-Sunyaev type vis- 
cosity prescription which was proposed for a Keplerian disc. 



And finally, the entropy generation equation : 



ds 



uT-- = Q^ -Q-, 
dx 



(Id) 



where, s is the specific entropy of the flow, T is the local temper ature. Q"*" and 
Q~ are th e heat gained and lost by the flow , and are given by (jChakrabartil . 



19961 : [Dad, l2007l : LShapiro fc Teukolskvl . Il983l ) 



and 



Q- 



--x[ga +7M ) — 



l3Sia^ 



ux 



3/2 1 



X 



Here g = In+i/In, n = 1/(7 - 1), In = (2"n!)V(2n + 1)! flMatsumoto et al. 



19841 ). and 7(= 4/3) is the adiabatic index. Presently, we consider only syn- 
chrotron cooling. In the above equation, /3 is the cooling parameter, and Si is 
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the synchrotron coohng term which is independent of the flow variables and 
is given by, 

where, e is the electron charge, rrie is electron mass, is the accretion rate 
in units of Eddington rate, Mq is solar mass, and for fully ionized plasma 
/i = 0.5. The suffix H = =f' represents the quantities in the pre/post shock 
disc region. It is to be borne in mind that in absence of shock m_|_ = m_, 
therefore 5"+ = S^. Due to the uncertainties of the realistic magnetic field 
structure in the accretion disc, we have assumed stochastic magnetic field. 
The ratio between the magnetic pressure and the gas pressure is represented 
by 1]. The magnetic field strength is estimated by assuming partial equipar- 
tition (?7<1) of the magnetic pressure with the gas pressure. In this paper, 
we have ig nored bremsstrahlung cooling, since i t is a very inefficient coo l- 



ing process (IChattopadhyay fc Chakrabartil. l2000l: iDas fc Chakrabartil . 120041 ). 



The expression for bremsstrahlung cooling (IRybicki fc Lightmanl . 119791 ) in ver- 
tical equilibrium is given by, 

Qs- 



ux^/'^{x — 1) ' 
where 

^ 2.016x10^^° 
' ~ Anml ^ 2A^'* 2GMqC 

where, rup is the proton mass and A;^ is the Boltzmann constant. For identical 
accretion rates 

S- 

= 3.26xl0^xr/. 

Therefore, it is quite evident that the synchrotron cooling is much stronger 
than bremsstrahlung. However, bremsstrahlung photons may interact with the 
accreting gas itself and in that sense bremsstrahlung may be important. Such 
complicated situation is not addressed in the present paper. We have also not 
considered inverse-Compton, since that will require a proper two temperature 
solution which is also beyond the scope of the present effort. 

In the present paper, we have chosen m_ = 0.1 and rj = 0.1 as the represen- 
tative case, until stated otherwise. 

Under the adiabatic assumption for the jet, the momentum balance equation 
can be represented in the following integrated form: 

^^■ = r? + ""? + ^-2(;^' 

where, Sj and Xj are the specific energy and angular momentum of the jet, 
respectively. Other fiow variables are the jet velocity (vj) and sound speed 
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(aj). Furthermore, Xj[= (xcb + xfw)/'^] and rj[= {x'j +1/'^^)^/'^] are the cyhn- 
drical and spherical radius of the je t streamhne. The funct i onal form of the 
coordinates of CB and FW are [see, IChattopadhyay fc ~Ds2 fl2007h ]. 



X 



FW 



XCB 



2X^.rcB{rcB - I] 



1/4 



2) + J(A| - 2)2 - 4(1 - y^^) 



^ 2 
where, xcb and xpw are measured at the same height of jet streamline and is 

given by ycB = \/{r'^B - ^cb)- 

The integrated form of mass-flux conservation equation for the jet is given by, 

Mout = PjVjA, (26) 

where, Mout is jet outflow rate and pj is the local density of the jet. The jet 
cross-sectional area is given by, A = 2tc{xqq — xj^^). 



3 Accretion-Ejection solution 



It is well known that matter falling onto black holes have to cross one or 
more critical points depending on the absence or presence of shock transitio n 
dChakrabartil . Il996l : IChakrabarti fc Dasl . I2004J : IChattopadhvav fc Dad . l2007h . 
If the flow parameters allow shock transition then matter must cross the sonic 
horizon twice, once before the shock and then after the shock. The location 
of the latter is called the inner critical point {xd) and the former is known 
as outer critical point (xco)- In absence of dissipation, the energy {£) and 
angular momentum (A) of the flow is conserved, and therefore Xd and/or 
uniquely obtained in terms of £ and A, and consequently all possible 
flow solutions. S and A do not remain conserved along a dissipative flow and 
therefore critical points cannot be determined uniquely. To obtain solutions 
of a dissipative accretion flow in a simpler way, one needs to know at least 
one set of critical point parameters {e.g., Xc, Ac). Fortunately, the range of 
(xri,\ri)s varies from (2r„ ^ x^ ~ 4r„,1.5 ^ Aw ^ A„,o), where Aw, A,,,, are 



the angular momentur n at the inner critical point and the marginally stab le 
orbit, respectively [e.g.. IChakrabartil fll989l . Il996h : IChakrabarti fc D ai tOOi \] . 
Here Vg is the Schwarzschild radius. Therefore for a viscous flow, it is easier to 
consider Xd and A^ as par ameters for solving t he flow equations, along with 
the v iscosity parameter a (IChakrabarti fc Dad . |2004 IChattopadhyay &: Dasl . 
20071 ). In presence of cooling, one should also supply the accretion rate at 
Xci in addition to {xd, Xd, C()- Presently, we fix accretion rate and vary /? to 
study the effect of cooling. Hence the existence of be obtained only in 

presence of a shock. 



7 



In this paper, we consider infinitesimally thin adiabatic shock, generally ex- 
pressed by the continuity of energy flux, mass flux and momentum flux across 
the shock, and is generally called Rankirie-Hug o niot (RH) shock conditions. 
Numerical simulations [e. a. . lEggum et al. I (ll985l ): lMolteni et al. I ( 1l994j . ll996al )] 
have shown that thermally driven outflows could originate from the hot inner 
part of the disc. When rotating matter accretes towards the black hole, cen- 
trifugal force acts as a barrier, inducing the formation of shock. At the shock, 
flow temperature rises sharply as the kinetic energy of the flow is converted 
into the thermal energy. This excess thermal energy may drive a signiflcant 
fraction of accreted material as outflows. Thus bulk properties such as excess 
thermal driving along z direction is a legitimate process for mass ejections. 



The modifled Rankine-Hugon iot shock conditions in presence of mass loss are 



Chattopadhyay fc Dad (120071 ) . and references therein] 



out 



MJI 



Rrn) 1 



n. = n_ 



(3) 



Assuming the jet to be launched with the same speciflc energy, angular mo- 
mentum and density as the post-shock disc, the expr ession for relative mass 
outflow rate is given by (IChattopadhvav fc Dad . l2007h . 



Mout/M- 



RVj{Xs)A{Xs 



An 



2 ^^/^ / ^ 



where, the compression ratio is deflned as i? = Since, the information 

of Rrn is in the shock condition itself, we need to solve accretion-ejection 
equations simultaneously. The method to do so is as follows: 

(a) we assume Rm = 0, (m_ = rh+), and with the supplied values of {xd, 
Xci, a, P) we integrate Eqs. (la-d) outwards along the sub-sonic branch of 
the post-shock region. Equation (3) is used to compute the pre-shock flow 
quantities, which are employed to integrate outwards to flnd the location of 
Xco- The location of the jump for which Xco exists is the virtual shock location 

«)■ 

(b) Once x'^ is found out, we assign £j = £{x'g) and Xj = \{x'J to solve the jet 
equations and compute the corresponding i?^. 

(c) We use this value of Rm in Eq. (3) and again calculate the shock location. 

(d) When the shock locations converge we have the actual shock location {xs), 
and the corresponding R^ is the mass outflow rate. 

In other words, we are launching jets with same £, A, and p as that of the 
shock. 



Presently, we consider viscosity and synchrotron cooling process as the source 
of dissipation in the flow. Viscosity reduces the angular momentum, while 
increases the energy as the flow accretes towards the central object. Cooling 
process on the other hand, decreases the flow energy inwards while leaving 
the angular momentum distribution unaffected. For proper understanding of 
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Fig. 2. £{x) with x is plotted for /? = (dashed), 0.01 (dotted) and 0.036 (soUd). 
Other parameters are (<Scj, Acj)=(0. 00182, 1.73) and a = 0.001. 

the effect of viscosity and coohng on determining mass outffow rates we need 
to fix {S, A) at some length-scale (around inner or outer boundary), and then 
vary a and /3. 

As Xci is very close to the horizon, fixing {Sd, Xd) at Xd is almost equivalent 
to fixing the inner boundary flow quantities. In Fig. 2, we plot £{x) with x for 
P — (dashed), 0.001 (dotted) and 0.0036 (solid), where the inner boundary 
flow quantities are Aci)=(0. 00182, 1.73) and a — 0.001. For the coohng free 
solution (dashed), the energy of the flow increases inwards due to viscosity. For 
solutions with significant cooling (dotted, solid), the increase in energy due to 
viscous heating is completely over shadowed, causing the energy to decrease 
towards the black hole. Increase in cooling efficiency signifies, matter with 
higher energies at the outer boundary, falls into the black hole with identical 
Sci- If standing shocks form, then under these circumstances energy at the 
shock will increase with (3. In the following, we discuss the role of viscous 
heating and synchrotron cooling in determining the mass outflow rate. 

In Fig. 3, we present a global inflow-outflow solution. In the top panel, the 
Mach number M of the accretion flow is plotted with log{x). The solid curve 
represents shock induced accretion solution. The inflow parameters are Xd = 
2.444, Xi = 1.75, a = 0.005, and /3 = 0.01 (for these parameters £d = 0.0018). 
In the lower panel, the outflow Mach number Mj is plotted with log{xj). In 
presence of mass loss, the shock forms at Xg = 21.64 denoted by the vertical 
line in the top panel, and the outflow is launched with energy and angular 
momentum at the shock {Ss,Xs = 0.00175, 1.766). The outflow is plotted up 
to its sonic point {xjc = 68.83), and the corresponding relative mass outflow 
rate is = 0.0816. 
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Fig. 3. Upper panel: Inflow Mach number (M = u/a) with log{x). The inflow 
parameters are Xd = 2.444, Aj = 1.75, a = 0.005, and /3 = 0.01 where, Xg = 21.64, 
Ss = 0.00175, As = 1.766 Xco = 166.57, A^ = 1.799. The dotted curve is the shock 
free solution. Lower panel: Outflow Mach number (Mj = vj/aj) with log{xj), the 
outflow critical point xjc = 68.63 {rjc = 270.8), and the jet coordinates at the base 
is given by xjb = 12.2 (r^b = 21.24). The relative mass loss rate is Rm = 0.0816. 




0.01 0.02 0.03 0.04 



Fig. 4. Variation of Rm with /3 for a = — 0.02 (left to right with da = 0.005). 
£■ = 0.0018 and A« = 1.75. 
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Fig. 5. Rrh is plotted with /3 for = -0.001^0.003 (right to left, dS^i = 0.001). 
Other parameters are Xd = 1.73 and a = 0.001. 



To present the global solution, Fig. 3 was obtained only for a set of input 
parameters, namely {Sd, Xd, a, P)- We would now proceed to find the explicit 
dependence of Rm on these parameters. In Fig. 4, we plot the mass outflow 
rates (Rm) with the cooling parameter (3, for a = — 0.02 (left to right for 
da = 0.005). All the curves are drawn for £d = 0.0018 and Xd = 1.75. Figure 4 
confirms o ur earlier investigation that R m decreases with increasing viscosity 



parameter (IChattopadhyay fc Dad . 120071 ). However, it may be noticed that for 



fixed a, Rm increases with /3. For a given a, the energy at the shock increases 
with P {e.g., Fig. 2), and since the post-shock region {i.e., CENBOL) is the 
base of the jet, the jets are launched with higher driving force. This causes 
Rm to increase with (3. It is to be noted, the two extreme curves {i.e., for 
a = 0.015, 0.02) on the right show that, for f5 = there is no outflow, but 
in presence of sufficient cooling steady jets reappear. As a is increased, Rm 
decreases due to the gradual reduction of sufficient dri ying at the jet base, and 



beyon d a critical a (say, acri) outflow rate vanishes (IChattopadhyay &: Das 



20071 ). For flows with a > the required jet driving could be generated by 
considering sufficiently high /3. In other words, to get steady outflows in the 
realm a > Ocri, there is a non-zero minimum value of (3 (say, (3m) corresponding 
to each a. Furthermore, for each a there is a cut-off in Rm at the higher end 
of P (say. Peri), since standing shock conditions are not satisfied there. Non- 
steady shocks may still form in those regions, and the investigation of such 
phenomena will be reported elsewhere. 

In Fig. 5, Rm is plotted with p for Sd = -0.001 (sohd), 0.0 (dotted), 0.001 
(big dashed), 0.002 (small dashed) and 0.003 (dash-dotted). Other parameters 
are Xd = 1-73 and a = 0.001. For a given P, mass outflow rate increases with 
£d- Higher £d corresponds to more energetic flow, and if these flows produce 
shock, we get higher Rm- On the other hand, even for same Sd, higher shock 
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Fig. 6. (a) Variation of Rrn with (i for Ac, = 1.73 (dotted) 1.75 (dashed) and 1.77 
(solid). £ci = 0.0018 and a = 0.001. (b) Variation of Rm, with Eg, for parameters 
same as Fig. 6a. 

energy is ensured with the increase of /9, and consequently higher Rfn are 
produced. The solutions corresponding to £ci = (dotted) and Sd = —0.001 
(solid) show that Rm^O as /3— >0. In other words, in presence of cooling, flows 
with bound energies at Xd may also produce outflows. Thus it is clear that 
shock energy plays an important role in determining the rate of mass loss from 
the disc. Previous studies of computation of mass outflow rates from inviscid 
and viscous disc showed that the angular momentum at the shock dictates 
the mass outflow rates, because higher angular momentum produces higher 
centrifugal driving for the jet. This lead us to investigate the role of angular 
momentum of the disc in determining the mass outflow rates, when cooling is 
present. 

In Fig. 6a, i?^ is plotted with /? for = 1.73 (dotted), 1.75 (dashed) and 1.77 
(solid), where — 0.00182, and a — 0.001 are kept flxed for all the curves. For 
negligible cooling {f3 ~ 0), higher angular momentum flow generates higher 
Rrn- As the centrifugal pressure produces the shock, which in turn drives 
the jet, it is not surprising that flows with larger angular momentum will 
produce higher i?^. Similar trend is maintained for nonzero f3. For a given X^, 
the energy at the shock (Sg) increases with /3. Thus the combined effects of 
centrifugal and thermal driving increase the mass outflow rate. We do see that 
there is a cut-off in R,n corresponding to each angular momentum at /3>/3cri. 
For lower angular momentum flow jScri is higher. To illustrate the effects of 
thermal driving and centrifugal driving of the jet, in Fig. 6b, we have plotted 
Rjh with Eg for Xci — 1.73 (dotted), 1.75 (dashed) and 1.77 (solid), for the same 
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Fig. 7. Rm is plotted with Sd for = 1.73 (dotted), Xd = 1-74 (dashed), Xd = 1.75 
(sohd). Other parameters are a = 0.001, and /? = 0.06. 

set of £ci and a as in the previous figure. It is to be remembered that Sg is not a 
new parameter but is calculated at the shock for the same range of (3 variation 
as in Fig. 6a. In the shaded region, Rm is higher for higher Xd- As long as the 
shock energy is similar, higher angular momentum results in greater centrifugal 
driving for the outflowing matter. However, lower ang ular momentum flo w can 



sustain higher energies across the shock [e.g., Fig. 3 of iDas et al. I (l2001al )]. For 
high enough Eg, the thermal driving starts to dominate over the centrifugal 
pressure, and results in higher Rm even for lower angular momentum flow. 

In Fig. 7, Rm is plotted as a function of £d, for various values of Xd = 1-73 
(dotted), Xd = 1-74 (dashed), Xd = 1-75 (solid). The other flow parameters 
are a = 0.001 and (3 = 0.06. This flgure distinctly shows that even if the ac- 
creting flow starts with unbound energy and produces shock induced outflow, 
signiflcant cooling closer to the black hole turns the unbound energy to bound 
energy. 



4 Astrophysical application 



In our solution procedure, we have employed three different constant parame- 
ters (3, 7] and rh to determine the cooling process. A cooling mechanism might 
depend on various other physical processes apart from its usual dependence on 
the flow variables. In general, rh regulates cooling, however to obtain a cool- 
ing free solution one needs to consider m = 0, which is meaningless. We have 
simplifled all such complications by introducing /? as a control-parameter for 
cooling. A simple inspection of Eq. (Id), shows that for a given set of (m, a, x), 
identical cooling rates may be obtained by rearranging the values of j3, rj and 
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m. It must be noted that, introduction of (3 and t] do not increase the param- 
eters of our solution, instead these are used to control the cooling efficiency 
and the magnetic field strength, about which there is no prior knowledge. In 
the previous section, we have fixed the values of m_ and t], and controlled the 
cooling term by jS.ln this section, we have fixed the value of (3 to unity, and 
allowed physical parameters, such as m_|_ and rj to dictate the cooling term. 



It is a matter of interest to estimate how much matter, energy and angular 
momentum enter into the black hole. In the present paper, the amount of mass 
fed to the disc is given by m_ . The rate at which matter is being accreted into 
the black hole and the rate of mass lo ss are self-consistently c o mput ed as m+ 
and (m_ — m+). It has been shown in IChattopadhyay fc Dad (120071 ) that the 
specific angular momentum of the fiow close to the horizon, is almost same 
as Xci- The actual value of S close to the black hole should be slightly higher 
than Sci- One has to quote the actual value of S close to the horizon. However, 
these numbers are obtained using pseudo-Newtonian potential and may not 
be consistent as general relativistic effects are important at such distances. 



We have applied our formalism to calculate the mass outfiow rates from two 
black hole candidates M87 and Sgr A*. M87 is supposed to harbour a super 
massive black hole [Mbh = 3xlO^ M^ (iFord et a/ I Il994h ]. The estimated 
accretion rate is M_~O.13M0yr~^ (IReynolds et al. I . Il996l ). The mass of the 
central black hole and the accretion rate of Sgr A* are Mrh = 2.6xl O^M(7) 



fISchodel et al. I . |2002[ ) and M_~8.8x lO-^Moyr"! flYuan et al I . l2002f ). The 



accretion disc around the black h ole in Sgr A* i s supposed to be radiatively 



in-efficient and of higher viscosity (IFalcke I . Il999l ). For both the cases we have 
set = 1, so the cooling mechanism is purely dictated by m and rj. To simplify 
further, we have chosen r] = 0.01 for both the objects. The accretion rates (in 
terms of Eddington rate) for M87 is given by m_ = 1.89 x 10~^ and that for Sgr 
A* is m_ = 1.47x10"^, therefore Sgr A* is dimmer than M87. With proper 
choice of a and Xd, and (see Table 1), we compute Rm (consequently 
m+) for both the objects mentioned above. The typical size of such a sub- 
Keplerian disc should be around a thousand Schwarzschild radii across the 
central object. Accordingly we have set the outer boundary at Xt = SOOr^,, 
and have provided the typical value of angular momentum at such distance 
(Aj-) for both the objects. For M87, the computed values of mass outfiow rate 
and shock location are Rm = 0.073 and Xs = 40.57. In case of Sgr A*, the 
estimated values of mass outfiow rate and shock location are Rm = 0.1049 
and = 14.415. 



Assuming the jet's luminosity is significant only at the lobes (where, the jet 
energy is mostly dissipated), the maximum luminosities of M87 and Sgr A* 
jets, estimated from the computed values of respective Rm, are given in Table 
1. Considering 10% radiative efficiency at the jet lobe the je t-luminosities for 



both M87 and Sgr A*, agree well with the observed values (IReynolds et al. 
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Table 1: Predicted values of Rm and jet power for M87 and Sgr A* 



Object 




M_ 


a 










Aj' 


Rfn 






Mq 


Mo/yr 




y 


cr„ 


MEdd 


y 


y 


% 


erg/s 


M87 


3.0 

xIQS 


0.13 


0.010 


2.367 


1.78 


1.75 

xl0~2 


40.57 


2.01 


7.3 


5.36 

xlO^^ 


Sgr A* 


2.6 

xlO^ 


8.80 

xlO"'^ 


0.015 


2.548 


1.71 


1.32 

xlO-4 


14.42 


2.44 


10.5 


5.2 

xlO^^ 



19961 : iFalcke fc Biermann I. Il999l). Mo r eover, the si ze of the comp u ted je t base 



for M87 is ~2x,~80rg. lJunor et al. I fll999f ) and iBiretta al. I fl2002t ) have 



estimated the base of jet to be less than lOOr^ from the central black hole, 
and probably greater than SOr^. Evidently our estimate of the jet base agrees 
quite well with the observations. There is no stringent upper limit of the jet 
base for Sgr A*, however, our computation gives a result which is acceptable 



in the literature (iFalcke 1 . 119991 ). We have also provided an estimate of angular 
mome ntum at Xr- Fo r Sgr A *, our estimated At is comparable with the re- 
sult of ICoker &: Melial (119971 ). However, no reliable estimate of At for M87 is 
currently available. In terms of physical units, various flow variables for M87 
are given by, Mout~0.009Moyr-\ M+~O.119M0yr-\ = S.lxlO^Vg g-\ 
a;s~3.61xl0^^cm, Acj~4.75xl0^^cm^s~^, and At~5.36x lO^^cm^s"^. Similarly 
for Sgr A*, Mout~9.1x IQ-^Moyr-i, M+~7.77xlO-^MQyr-i, = AxlO^^eig 



Xa 



^.llxlO^^cm, Ac»~3.96xl022cm2s-\ and AT~5.64xl022cm2s- 



In this paper, only sub-Keplerian matter dis t ributi on is chosen for the accre- 
tion d isc. However. IChakrabarti fc TitarchukI (119951 ) and lChakrabarti fc Mandal 
(120061 ) have shown that if a mixture of Keplerian and sub-Keplerian matter 
is chosen, then the spectral properties of the disc is better unde rstood. These 



assertions have been ratified for several black hole candidates ( ISmith et al. 



2OOII . I2OO2I ). Since matter close to the black hole must be sub-Keplerian, there- 
fore regardless of their origin, Keplerian and sub-Keplerian matter mixes to 
produce sub-Keplerian flow before falling onto the black hole. Such transition 
from two component t o single component flow has been shown by various au- 



ipc 

thors [e.g., Fig. 4b, of Das et al. (j2001bl )]. The region where such transition 



occurs may be called 'transition radius' (Xt). It must be noted that, Xt is 
treated as the 'outer edge' of the disc in our formalism described so far. The 
energy (St) and angular momentum (At) at X^ can then easily be expressed 
in terms of the accretion rate of the Keplerian co mponent {Mk) and the sub- 



Keplerian component [Msk) (IDas et al. 
and the net accretion rate being M 



2001bD 



Once Xt, St, Xt is known 
Msk + Mk, it is easy to calculate 
Rrn following our formalism. Thus, it is p ossible to predict Rm from t he sp ec- 
trum of the accretion disc, if formalism of IChakrabarti fc TitarchukI (119951 ) is 
applied on our solutions. 



15 



5 Concluding Remarks 



The main goal of this paper was to study how dissipative process e s affe ct the 



jet generation in an advective disc model. IChattopadhyay &: Dad (120071 ) have 
shown that mass outflow rates decrease with increasing viscosity parameter. 
In the present paper, we have investigated how the mass outflow rate responds 
to the synchrotron cooling. The general method of the solution (succinctly de- 
scribed in Section 3.) is to supply Xd, Xd, a, [3 and then integrate outwards to 
find the shock location (and consequently the mass outflow rate). Needless to 
say, once the above four parameters are fixed, the solution determines flow with 
unique outer boundary {i.e., at Xt). Of the four parameters, if a is increased, 
the solution corresponds to flow with higher angular momentum and lower 
energy at the outer boundary. On the contrary, when j3 is increased then the 
solution corresponds to higher energy but identical angular momentum flow 
at the outer boundary. Consequently, more energetic flows are allowed to pass 
through standing shock for higher /3, and hence stronger jets are produced. 
We have also shown that, if cooling efficiency is increased, then it is possible 
to produce jets even for those a-s for which i?^ is zero {e.g., Fig. 3). Fur- 
thermore, it has been shown that the jets are primarily centrifugal pressure 
driven even in presence of cooling. We notice that standing shocks in higher 
angular momentum flow do not exist for higher cooling efficiency, therefore 
steady jets are not produced. However, for higher /?, low angular momentum 
flow can generate high enough relative mass outflow rates. 

We have applied our formalism on a couple of black hole candidates, namely, 
Sgr A* and M87. Using the available accretion parameters of the above two 
objects as inputs, we have shown that one can predict observational estimates 
of jet power. More over, the typic a l size of the jet base (~ 2xs) also agrees well 



with observations. iLe &: Beckerl (120051 ) had dealt with these two particular 



objects, with their methodology which also involve shocked accretion disc. 
The methodologies of the present paper and the work of Le & Becker (2005) 
is quite different in the sense that, Le & Becker (2005) dealt with isothermal 
shock while our model is based on the adiabatic shock scenario. In Le & Becker 
(2005), the focus was on calculating the number densities and energy densities 
around an isothermal shock of an hot tenuous adiabatic rotating flow, by first 
order Fermi acceleration process. The energy lost at the isothermal shock, 
drives a small fraction of in falling gas to relativistic energies. With the given 
observational estimates of black hole mass, accretion rate etc of M87 and Sgr 
A*, they estimated the Lorentz factors of the jet. We on the other hand, have 
computed the thermally driven outflows from the post-shock disc, where the 
jets are launched with the local values (specific energy, angular momentum and 
density) of the disc fluid at the shock. With input values of black hole mass, 
accretion rate, and proper choice of viscosity parameter, inner sonic point etc 
we predict the shock location, the mass outflow rate. We check whether the 
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predicted values are within the accepted hmits or not. We do not estimate 
the terminal bulk Lorentz factor, since we believe one has to recast the whole 
framework into the relativistic domain as well as employ other accelerating 
processes {e.g., magnetic fields etc). One may wonder at the veracity of the 
two different processes employed to explain the observational estimates of jet 
quantities of M87 and Sgr A*, in other words, whether the jets are generated 
by post-shock thermal driving (we have not investigated magneto-thermal 
driving since this is only hydrodynamic investigation), or the jets are launched 
by particle acceleration processes. If one can observationally estimate the rate 
at which mass being ejected from the accretion disc, probably then one can 
ascertain the dominant effect behind jet generation. If it can be established 
that indeed the rate of mass los s is negligible compa red to the accretion rate 



then probably the formalism of iLe fc Beckerl (120051 ) is the more realistic jet 



generation mechanism. However, suffice is to say, various numerical simulation 
results do show (for non-dissipative as well as dissipative flows) that post-shock 
flow thermally drive bipolar outflows, and our effort has been to investigate 
how dissipative processes affect the relative mass outflow rates. 

In this paper we have only discuss ed formation of steady jets, since we have 



considered only stationary shocks. iMolteni et al. I (jl996bl ) have shown that 



the periodic breathing of the CENBOL starts when the post shock in-fall 
timescale matches with the Bremsstrahlung cooling timescale. Presently, we 
have considered dissipative processes which are more effective in determin- 
ing shock properties compared to Bremsstrahlung. Therefore, the dissipative 
processes considered in this paper, may trigger comparable or different shock- 



instab ilities in the disc than that has been reported earlier (iMolteni et al. 



1996bl ). Since, the jet formation is primarily controlled by the properties of 
the shock, any non-steady behaviour of the shock will leave its signature on 
the jet. In particular, a significant oscillation of the shock (both in terms of 
the oscillation frequency and its amplitude) may produce periodic ejections. 
We are studying dynamical behaviour of the shock in presence of viscosity and 
synchrotron cooling using fully time dependent simulation and results will be 
reported elsewhere. 
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